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The independent eigenstates of the total orbital angular momentum 
operators for a three-body system in an arbitrary D-dimensional space are 
presented by the method of group theory. The Schrodinger equation is re- 
duced to the generalized radial equations satisfied by the generalized radial 
functions with a given total orbital angular momentum denoted by a Young 
diagram [/i, z^, 0, . . . , 0] for the SO{D) group. Only three internal variables 
are involved in the functions and equations. The number of both the func- 
tions and the equations for the given angular momentum is finite and equal 
to — V + 1). 

I. INTRODUCTION 

From the very early stage of the progress in quantum mechanics in the real three- 
dimensional world, it has been pointed out that the essence of these theories would be 
easily understandable if their mathematics is constructed in the non-relativistic hyper- 
space worlds The mathematical tools for generalization of the orbital angular mo- 
mentum in an arbitrary D-dimensional space have been presented Recently, the 
D-dimensional Coulombic and the harmonic oscillator problems in a two-body system 
have been studied in some detail by many authors 

Exact solutions played very important roles in the development of physics. The exact 
solutions of the Schrodinger equation in the real three dimensional space for a hydrogen 
atom and for a harmonic oscillator were important technical achievements in quantum 
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mechanics which provided strong evidence in favor of the theory being correct, at 
least as far as atomic physics is concerned. The next simplest atom is the helium atom, 
for which the Schrodinger equation cannot be solved analytically, but only numerically 
2^^8[. In the numerical calculation, one of the main difficulties is how to separate the 



global rotational degrees of freedom. 

In our previous paper we separated completely the global rotational degrees of 
freedom in the Schrodinger equation for an A^-body system in the real three-dimensional 
space from the internal ones. We have determined a complete set of {21 + 1) independent 
base functions for a given total orbital angular momentum /, which are the homogeneous 
polynomials in the components of coordinate vectors and do not contain the Euler angles 
explicitly. Any function with the given angular momentum / in the system can be ex- 
panded with respect to the base functions, where the coefficients are the functions of the 
internal variables, called the generalized radial functions. The generalized radial equa- 
tions satisfied by the functions are established explicitly [^. For the typical three-body 



system in the real three dimensional space [pO| , pT| , such as a helium atom [0,0 and 
a positronium negative ion |34], the generalized radial equations [3^ have been solved 
numerically with high precision. 

With the interest of higher dimensional field theory recently, we have a try to generalize 
the study of the D-dimensional two-body system to the D-dimensional three-body system. 
The purpose of this paper is, for a three-body system in an arbitrary D-dimensional space, 
to find a complete set of independent base functions with any given total orbital angular 
momentum and to reduce the Schrodinger equation with a spherically symmetric potential 
V to the generalized radial equations, where only three internal variables are involved. 
Any function with the given angular momentum in the system can be expanded with 
respect to the base functions. It provides a possibility to calculate numerically the energy 
levels of the three-body system in Z)-dimensions with high precision. 

From the viewpoint of mathematics, the separation of the global rotational degrees of 
freedom from the internal ones is a typical application of group theory to physics. The 
properties of the independent base functions for a given total orbital angular momentum 
would be more clear if they are constructed in arbitrary D-dimensional space than that in 
the real three dimensional space. The total orbital angular momentum for a three-body 



2 



system in a D-dimensional space is described by an irreducible representation denoted by a 
Young diagram with one or two rows. For the real three-dimensional space, the rotational 
symmetry group is SO (3) group, and its only irreducible representations denoted by the 
Young diagrams with two rows are [/, 1], which are equivalent to the representations 
denoted by the one-row Young diagrams [/,0], respectively. This is the reason why the 
angular momentum can be described by only one quantum number I for the real three- 
dimensional space. 

This paper is organized as follows. After separating the motion of the center of mass 
by the Jacobi coordinate vectors in Sec. 2, we review in Sec. 3 the generalization of the 
orbital angular momentum operators and the properties of the spherical harmonics 
and the harmonic polynomials [[^] for a two-body system in D dimensions. In Sec. 4 we 
will define the generalized harmonic polynomials for a three-body system in D dimensions 
and prove that they constitute a complete set of independent base functions for a given 
total orbital angular momentum in the system. The generalized radial functions are 
defined and the generalized radial equations are derived in Sec. 5. Some conclusions will 
be given in Sec. 6. 



For a quantum A^-body system in an arbitrary D-dimensional space, we denote the 
position vectors and the masses of particles by and by m^, k = 1,2, . . . , N, respec- 
tively. M = J2k '^k is the total mass. The Schrodinger equation for the A^-body system 
with a pair potential V , depending upon the distance of each pair of particles, |rj — r^l, is 



where v?j. is the Laplace operator with respect to the position vector r^. For simplicity, 
the natural units Ti = c = 1 are employed throughout this paper. The total orbital angular 
momentum operators Lab in D dimensions are defined as 



II. SCHRODINGER EQUATION IN D DIMENSIONS 



1 ^ 



(1) 



k=l 




a,b = 1, 2, . . . 



(2) 



where r^a denotes the ath component of the position vector r/j. 
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Now, we replace the position vectors by the Jacobi coordinate vectors Rj-: 

1 < J < (iV - 1), M, = mfc, Mn = M, (3) 

fc=i 

where Rq describes the position of the center of mass, Ri describes the mass-weighted 
separation from the second particle to the first particle. R2 describes the mass-weighted 
separation from the third particle to the center of mass of the first two particles, and 
so on. An additional factor v^M is included in Rj for convenience. The mass-weighted 
factors in front of the formulas for Rj are determined by the condition 

N N-1 



k=l j=0 

One may determine the factors one by one from the following schemes. In the center- 
of-mass frame, if the first j particles coincide with each other and the last (A^ — j — 1) 
particles are located at the origin, the factor in front of Hj is determined by 

n = T2 = ■ ■ ■ = Tj = -mj+iTj+i/Mj, ^k^l = R-j- (4) 

k=l 

A straightforward calculation by replacement of variables shows that the Laplace op- 
erator in Eq. (1) and the total orbital angular momentum operator Lab in Eq- (2) are 
directly expressed in R^: 

N N-1 
V' = E "^fc Vr, = E Vr,, 

N-1 r n n ^ \'~') 



N-1 ( 
Lab = ~^ E 1 
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dRjb dRja J 

In the center-of-mass frame, Rq = 0. The Laplace operator (5) obviously has the 
symmetry of the 0{ND — D) group with respect to (A^ — 1)D components of (A^ — 1) 
Jacobi coordinate vectors. The 0{ND — D) group contains a subgroup SO(D) xO(A^ — 1), 
where SO{D) is the rotation group in the D-dimensional space. The space inversion and 
the different definitions for the Jacobi coordinate vectors in the so-called Jacobi tree 
can be obtained by 0(A^ — 1) transformations. For the system of identical particles, the 
permutation group among particles is also a subgroup of the 0(A^ — 1) group [p3]. 



It is easy to obtain the inverse transformation of Eq. (3): 



M. 



mjMj_i 



-[ 1/2 N-1 
k=j 

1/2 j-2 

+ E 

i=k 



MkMk+i 



ni/2 

±li — 



R 



(6) 



Thus, the potential is a function of R,- • Rjt and is rotationally invariant. 



III. HARMONIC POLYNOMIALS IN D DIMENSIONS 

In the center-of-mass frame, Rq = 0. Hence, for a two-body system there is only one 
Jacobi coordinate vector Ri, which will be denoted by x for simplicity: 
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Louck 1^,^ introduced the hyperspherical coordinates 

Xi = r cos 6*1 sin 6*2 . . . sm9£)_i, 

Xo = r sin 6*1 sin 60 ■ ■ ■ sin 6'n_i , 

' (8) 
Xk = r cos6'fc„i sin 6'^ . . . sin 60-1, 3 < k < D — 1, 

xu = rcosOo^i- 

The spherical harmonics ... l-^ in D dimensions are the simultaneous eigenfunc- 
tions of the commutant operators L^: 

T2- t2- / 1 ^ .,^k-ln 9 LLi 1 

LlYl^_,,...,i,{ei . . . en-i) = hih + k - J^i • • • ^^-1)' (10) 

/ = Id-1 = 0, 1, . . . , /fc = 0, 1, . . . , lk+1, h = —h, —h + 1, • • • , ^2 — 1, ^2, 
where = Lf^.^, < r < 00, -vr < < tt, < 0^ < vr, and 2 < A; < D - 1. The 
volume element of the configuration space is BJT4 



D D~l 



\{dXj = r^-Hr\{ {sme^y-Uej. (11) 
i=i i=i 



Through a direct calculation by replacement of variables, one obtains 

= ^ Ar^-i A _ ^ (12) 

Due to the spherical symmetry, the wave function can be expressed as 

Jx) = • • • ^^-i)' (13) 

and the D-dimensional Schrodinger equation (1) for a two-body system in the center-of- 
mass frame reduces to the radial equation 

- ^ ^h iir) = -2[E- Vir)] 0,(r). (14) 

Bander and Itzykson introduced the harmonic polynomials in D dimensions 

3^L.,---,Jx) = ^X_.--a(^i---^^-i) (15) 

to avoid the angular functions Yl^_^ ... ^^(6*1 . . . 6d-i)- 3^/o_2 ■■■ 'i(-'^) ^ homogeneous poly- 
nomial of degree I in the components of x and satisfies the Laplace equation 

Vx3^L,,., Jx) = 0. (16) 

The number of linearly independent homogeneous polynomials of degree I in D com- 
ponents of X is A^(/) = {I + D — iy./l\{D — 1)!. The Laplace equation (16) gives 
A^(/ — 2) = {I + D — 3)1/ {I — 2)\{D — 1)\ constraints. Hence, the number of the harmonic 
polynomials ^/^^...^^(x) of degree / as well as the number of the spherical harmonics 
■■■Zi(x) D dimensions is 

m - - 2) ^ ^"^^-^^^!;r'^'' - ^.(M, .... 0]). (17) 

(i£)([/, 0, . . . , 0]) is the dimension of the irreducible representation of SO{D) denoted by the 
one-row Young diagram [/, 0, . . . , 0]. [/, 0, . . . , 0] describes the symmetric traceless tensor 
representation. In fact, any polynomial in the components of one vector x has to belong 
to a symmetric representation. 

Due to the spherical symmetry, one only needs to write the explicit form of the highest 
weight state 

yl..^ii^) = Niixi + tX2y. (18) 



where Ni denotes the normahzation factor. The partners of iV(...^/(x) can be simply gen- 
erated by rotation. Now, the solution to the Schrodinger equation in the center-of-mass 

frame can be re-expressed as 

= i?Kr)3^i...,,(x), (19) 
and the radial equation is easy to be derived: 

1 ^-r^-i|-it:Kr) + -^Ri{r) ^ -2[E - V{r)] Ri{r). (20) 



Recall Ri{r) — r~''(f)i{r). Eq. (20) coincides with Eq. (14) but the angle variables do not 
appear explicitly in calculation. 

The number (17) of the harmonic polynomials 3^/^_2 .../i(x) of degree / can be under- 
stood from another viewpoint. After removing those homogeneous polynomials in the 
form r^/(x), where /(x) is a homogeneous polynomial of degree {I — 2), Eq. (17) shows 
the number of the remaining linearly independent homogeneous polynomials of degree I 
in the components of x. Therefore, the harmonic polynomials 3^/o_2,-,;i(x) construct a 
complete set of linearly independent base functions for the homogeneous polynomials of 
degree / in the components of x, excluded those in the form of r^/(x). 

IV. THREE-BODY SYSTEM IN £)-DIMENSIONS 

For a three-body system, in the center-of-mass frame there are two Jacobi coordinate 
vectors Ri and R2, which will be denoted by x and y, respectively: 

1/2 



mim2 

TTli + 1712 



{r2 - ri} , y 



"1 1/2 

(mi + 7722)7723 [ TTliTi + 7?72r2 



mi -|- m2 + ms 



mi -|- m2 



= Vx + Vy, 

-ix) , j{y) .id d ] .{ d d ] (21) 



The Schrodinger equation (1) reduces to 

{ + v^} y) ^-2{E- y(6, 6, 6)} y), 
= X • X, ^2 = y ■ y, = X ■ y. 



(22) 
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where are the internal variables. Since Eq. (22) is rotational invariant, the total orbital 
angular momentum is conserved. The wave function ^(x, y) with the given total angular 

momentum has to belong to an irreducible representation of SO(i^). In the traditional 
method, one calculates the wave function by the Clebsch-Gordan coefficients: 

As usual, 3^/o_2, .■■,Zi (x) and 3^/^ ^^..-^vjc/) ^-^e called the partial angular momentum states, 
and their combination is called the total angular momentum state, which is a homogeneous 
polynomial of degrees I and V in the components of x and y, respectively. 

There are three problems. First, what kinds of representations (or total angular mo- 
mentum L) appear in the Clebsch-Gordan series for decomposition of the direct product of 
two representations denoted by one-row Young diagrams [/, 0, . . . , 0] and [/', 0, . . . , 0]? This 
problem has been solved in group theory by the Littlewood- Richardson rule and traceless 
conditions. A new character is that the representations denoted by two-row Young dia- 
grams appear in the Clebsch-Gordan series for a three-body system when D > 3. Those 
representations denoted by the Young diagrams with more than two rows could not appear 
because there arc only two Jacobi coordinate vectors. For simplicity we denote a one-row 
or two-row Young diagram by [/x, v] = [/x, i/, 0, . . . , 0]. Hence, we have the Clebsch-Gordan 
series: 

[Z,O]®[r,O]~00 [l + l' -s-2t,sl (24) 

where n is the minimum between / and /'. The representations with t = arc calculated 
by the Littlewood-Richardson rule, and the remaining are calculated by the traceless 
conditions. The dimension of a representation denoted by a two-row Young diagram is 

Mil^.^]) = (^ + 2/i-2)(D + /x + z/-3)(D + 2t/-4)(/^-t/ + l) 

{D + ^1 - A)\{D + u - h)\ (25) 
^ {ji + l)\u\{D - 2)\{D - 

When D = 4, the representation denoted by a two-row Young diagram reduces to a 

direct sum of a selfdual representation [{S)^,^] and an antiselfdual one [{A)^,v\. Their 

dimensions are equal to half of di{\pi,v\) given in Eq. (25). When D — due to 

the traceless condition, the only representations with the two-row Young diagrams are 

representations [/x, 1], which are equivalent to that with the one- row Young diagrams [/i., 0], 
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respectively. Eq. (25) still holds for D = 3. Second, how to calculate the Clebsch-Gordan 
coefficients? The calculation must be very complicated. We will avoid the difficulty by the 
method of determining the highest weight states directly. Third, how many base functions 
are independent for a given total orbital angular momentum such that any wave function 
with the same angular momentum can be expanded with respect to the base functions 
where the coefficients are the functions of the internal variables. We are going to solve 
the last two problems by group theory. 

Let us sketch some necessary knowledge of group theory. From the representation 
theory of Lie groups [^^-|38| , the Lie algebras of the S0(2n+1) group and the SO(2?2) group 
are and Dn, respectively. Their Chevalley bases with the subscript j, 1 < j < n — 1, 
are same: 

= L(2j-l){2j) - L(^2j+l){2j+2), 

= {L(2j)(2j+i) - iL(2j-i){2j+i) - iL(^2j){2j+2) " L(^2j~i)(2j+2)) /2, (26a) 

^3 = (^(2i)(2i+l) + «^(2j-l)(2j+l) + «^(2i)(2j+2) " ^(2i-l)(2i+2)) /2. 

But, the bases with the subscript n are different: 

Hn = 2L(2n-l)(2n); 

En = L(^2n){2n+1) — ^-^(2n-l)(2n+l) , {26b) 
Fn = L[2n)(2n+l) + ^-^(2n-l)(2n+l) , 

for S0(2n + 1), and 

Hn = -^(2n-3){2n-2) + -^(2n-l)(2n) , 

En = {L{2n-2){2n~l) " '^-^(2n-3)(2n-l) + ^-^^(2n-2)(2n) + -^(2n-3)(2n) ) /2, (266) 
En = {L{2n-2){2n-l) + ^-^(2n.-3)(2n-l) ~ iL(2n~2){2n) + -^(2n-3)(2n) ) /2, 

for S0(2n). span the Cartan subalgebra, and their eigenvalues are the components of 
a weight vector m = (mi, . . . , m„): 

iffc|m) = mfc|m), 1 <k <n. (27) 

If the eigenstates for a given weight m are degeneracy, this weight is called a multiple 
weight, otherwise a simple one. Ek are called the raising operators and Ek the lowering 
ones. For an irreducible representation denoted by a Young diagram [/xi, /i2, . . .] of SO(D), 
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fJ'j > there is a highest weight M = (Mi, M2, . . .), which must be simple: 
Mj = Hj - /ij+i, 1 < j < ^ - 2, 

Mn-l = lin-l - fJ'n, Mn = 2/i„, for S0(2n + 1), 

Mn-i = /in-1 — fJ^n, Mn = yUn-1 + yUn, for selfdual representation in S0(2?t,), 

Mn-i = /in-i + A^n, Mn = " /^n, for antiselfdual representation in S0(2n). 

(28) 

We are not interested here in the spinor representations where Mn is odd for S0(2n + 1) 
and Mn-i + Mn is odd for SO(2?2). For a given irreducible representation [/ii,/i2, . . .] of 
SO{D), we only need to consider the highest weight state |M), which satisfies 

Hk\M)=Mk\M), Ek\M) = 0, l<k<n, (29) 

because its partners can be calculated by the lowering operators Fk- In this paper the 
highest weight state will simply be called the wave functions with the given angular 
momentum [/i, u] for simplicity. 

Now, we return to our problems. Recalling the Clebsch-Gordan series in Eq. (24), we 
can rewrite Eq. (23) for the highest weight state M: 

yS'''\^,y) = i: 3^L(x)3^M-m(y)(^,m,/',(M-m) I [(/ + /'-,-2t),.],M), (30) 

m 

where the subscripts of the harmonic polynomials are changed to the weights for simplicity. 

''^'*(x, y) is the highest weight state of the representation [(/ + /' — s — 2t), s]. It is a 
homogeneous polynomial of degrees I and I' in the components of x and y, respectively. 
Generally speaking, some ''''*(x, y) may be expressed as a sum where each term is a 
product of an internal variable C,j and a homogeneous polynomial /(x, y) of lower degree 



(see p. 042108-5 in p9|]). Since 3^M'^'*(X'y) "^i^^ ^e used as a base function for the wave 
function with a given angular momentum and the combinative coefficient is the function 
of the internal variables, in this meaning, the base function in the form of ^j/(x, y) is not 
independent, and we should find out the independent and complete base functions for any 
given angular momentum [/i, z/]. In the following we are going to prove y^^'^'^ ''^''^'"(x, y) 
and their partners, where I = q, I' = fi + u — q, s = u, t = 0, and < q < constitute a 
complete set of independent base functions for the total orbital angular momentum [fi, v\ . 
In other words, those total angular momentum states 3^m '"^'^l^' y) "with t > are not 
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independent, where the sum of the partial angular momentum quantum number / and /' 
is larger than fi + u for the total angular momentum [/i, u] . 

The highest weight for the representation [/i, u] is M = (/i — i/, z/, 0, . . . , 0). Removing 
the normalization factor in 3^m^^'^~'^''''^'°(^) y)) which is irrelevant here, we can determine 
the explicit form for y^^'^'^~'^^''^'^{'x.,y) according to its orders in the components of x 
and y and the property of the highest weight state (29), and denote it by the generalized 
harmonic polynomial (^^'^(x, y): 



y^,+u-,),.,o^^^^^^ 0<z.<g</., (31) 



Xi2 = xi+ ix2, X34 = X3 + ix4, Y12 = yi + iy2, >34 = 1/3 + iy^- 

The formula for (x, y) holds for D = 3 (X4 = ^4 = 0, z/ = or 1) and D > A. 

When D = 4 we denote the highest weig ht states by Qf)^''(x,y) and Q(^)'^^(x,y) for the 
selfdual representations and the antiselfdual representations, respectively: 

^34 = ^3- ixi, Y34 = 1/3 - iyi- 

The generalized harmonic polynomial Q^'^(x, y) is a homogeneous polynomial of degrees 
q and (/i + z/ — g) in the components of x and y, respectively. It is a simultaneous 
eigenfunction of Vx; Vy^ Vx ■ Vyi the total angular momentum operator [see Eq. 



v^Qr(x, y) = VyQr(^' y) = Vx ■ VyQr(^' y) = 0' 

L2QM^(x,y)=C2([^,z/])gf(x,y), (33) 

C2{[fi, u]) = fi{fi + D~2) + u{u + D- 4), 
where C2{[fi, t^]) is the Casimir calculated by a general formula (see (1.131) in Ref. |P8|| ). 
The parity of (x,y) is obviously (-1)^+^ 

It is evident that (^^'^(x, y) do not contain a function of the internal variables as a 
factor, neither do their partners due to the rotational symmetry. Therefore, Q^^(x, y) are 
independent base functions for the given angular momentum described by [/i, u] . Now, 
we are going to prove that {fi — u + 1) base functions (^^'^(x, y) where v < q < jJ. are 
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complete for the angular momentum u] . That is, Q^^''"^^ (x, y) with Q <l — ji < q < n 
and their partners construct a complete set of linearly independent base functions for the 
homogeneous polynomials of degree / in the components of x and y, excluded those in 
the forms of ^j/(x, y), where /(x, y) is a homogeneous polynomial of degree (/ — 2). 

The number of linearly independent homogeneous polynomials of degree / in the com- 
ponents of X and y is 



Md{1) 



I 



2D-1 

After removing those polynomials in the form ^jf{^, y), the number M(/) reduces to K{1): 
Kd{1) = Md{1) - SMd{1 - 2) + 3Md{1 - 4) - Md{1 - 6) 

(i^2D-7V (34) 
= 4(Z + D - 3) [2l{l + 2D-6) + {D- 2){2D - 5)] _ , 

when I + 2D > 7, which only excludes one case of Z = and D — 3, where K3{0) — 1. 

On the other hand, the number of Q*^^^'^^ (x, y) with 0<l — ijL<q<ijL and their 
partners can be calculated directly from Eq. (25): 

J2 (2// - / + l)dD{[f^, {I -//)]) = Kd{1). (35) 

1/2<IJ.<1 

Eqs. (34) and (35) are checked by Mathematica. Thus, we have proved that {/j, — u + 1) 
base functions Q^''(x, y) where < u < q < /j, are independent and complete for the 
angular momentum [/x, u]. Any function with the angular momentum [n, v\ in the system 
can be expanded with respect to the base functions (^^'^(x, y), where the coefficients are 
functions of internal variables. 

Prom Eq. (30), for a given total orbital angular momentum [//, i^] there are infinite 
number of wave functions y-j combined from different partial angular 

momentum states. Now, what we have proved is that only a finite number of partial 
angular momentum states {t = 0) are involved in the complete set of independent base 
functions (5g^(x, y) for a given total orbital angular momentum [/x, v]. 

V. GENERALIZED RADIAL EQUATIONS 

In the preceding section we proved that any function with angular momentum [/i, z/] 
in the quantum three-body system of D dimensions can be expanded with respect to the 
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base functions (5^''(x, y) 

^te'^kx, y) = E V'r(^i' e3)gr(^' y)' (36) 

q=v 

where the coefficients Cs) ^'^^ called the generalized radial functions. When 

substituting Eq. (36) into the Schrodinger equation (22), the main calculation in the 
derivation is to apply the Laplace operator (21) to the function ^!^\yi^ y). The calcula- 
tion consists of three parts. The first is to apply the Laplace operator to the generalized 
radial functions '^^'^(Ci, ^2) Ca); which can be calculated by replacement of variables: 

V'<^(6, 6, 6) = [^iidl + Ai,dl + 2D (% + a^) + (6 + 6) dl ^^^^ 

+4e3(% + %)%}v^r(^i'^2,e3), 

where 9^ denotes d/d^ and so on. The second is to apply it to the generalized harmonic 
polynomials Qf^^ (x, y ) . This part is vanishing because Q^'' (x, y ) satisfies the Laplace 
equation. The last is the mixed application 

2 { (a^.^f ) 2x + ) y } ■ VxQr + 2 { ) 2y + ) x} ■ VyQ^- (38) 

Prom the definition (31) for (5^'^(x, y) we have 



(39) 



X • VxQf = gQf, y • VyQf = ()" + ^ - <i)Q7 

y ■ VxQr = (//-? + i)Q^-i, X • VyQr = (? - ^ + i)Q^+i. 

Hence, wc obtain the generalized radial equation, satisfied by the {/j, — + 1) functions 

V'r(^i,e2,e3): 

+ 4g%<'^ + 4(/. + - g)%^^^- + 2(/i - g)%<^i + 2(g - i/)%d ^^^^ 

= -2 (^-l^)V'f, 

where \j'^'4>q^ is given in Eq. (37). Only three invariant variables ^1, ^2 and ^3 are 
involved both in the equations and in the functions. When D = 4, Eq. (40) holds for the 
generalized radial functions either in [(<S')/x, u] or in [(A)/x, z/], because two representations 
incorporate to one irreducible representation of the 0(4) group when the space inversion is 
considered. When D = 3, the equations for the functions in [//, 0] and in 1] are different 
although two representations [/x, 0] and [/x, 1] are equivalent, because the functions have 
different parity. 
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At last, we discuss rotational variables and the volume element of the configuration 
space. We fix the body-fixed frame such that x is parallel with its Dth axis, and y is 
located in its {D — 1)D hyperplane with a non-negative (D — l)th component. That is, 
in the body-fixed frame, the nonvanishing components of two Jacobi coordinate vectors 
x' and y' are 

^ = e}/^ y'n-, = {^2-ej^iy^\ y'n = i.ii'"- (41) 

Let R = eSO(D) rotate the center-of- mass frame to the body-fixed frame: 

i?*^^^ = -Rl2(^l)-R3l(^2)-R43(^3)-R54(^4) ■ ■ ■ Rd{D-1){0 D-l) , 

= Rl2MRsi{^p2)Riz{V-i)RbA{V4) ■ ■ ■ R{D-1){D-2){VD-2), (42) 

i?x' = X, Ry' = y, 

where, for example, RuiO) is a rotation on the hyperplane with the first and the second 
axes through 9 angle: 



^ cos ^ — sin ^ ^ 



Ri2{e) 



V 



sin 9 cos 9 
1d-2 

{D — 1) 9j and {D — 2) ipk are the rotational variables, called the generalized Euler angles. 
Through a straightforward calculation, we obtain 



Xi + 1X2 = ^y^e^^i sin 9a-, Xs + ix4 = ^i '^ (cos 6*2 sin 9^ + i cos ^3) sin 9a, 



D-l 



1/2 



D-l 



a=2 



D-l 



a=1 



yi + iy2 = izil^'^^'^ n sin 9a + (6 - ^37^1) 

0=2 

( D-2 D-2 /a \ / D-2 

X \ iW sin ^ cos 9a+i cos ( H sin 6*6 j J| sin 



a=l 



a=l 



^6=2 
D-l 



ys 



+ m = ^3^1 ^^"^ (cos 612 sin 93 + i cos 6^3) ]J sin 9a + - Cl/^i 



a=4 



1/2 



{_D~2 D-2 
— COS (fii sin ^2 n ^i^ V"- + (^"^^ ^2 cos 9s — i sin ^3) cos 9?2 Yl ^in (/'a 
a=2 a=3 
D-2 /a \ / D-2 

-I- (cos ^2 sin ^3 -I- i COS ^3) ^ cos ^a+i cos j JJ sin ^6 sin(/?c 

a=3 \6=4 / \c=a+l 

where na=6+i Fa — ^- The volume element of the configuration space is 



(43) 



D ^ 

Yl dxjdyj = - (Ci6 - it) 



(D-3)/2 



diidi2disl[ (sin^,)'"'c^^,- n' {sm^k)"'' d^^- (44) 
j=i k=i 
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VI. CONCLUSIONS 



After separating the motion of center of mass, we have defined the homogeneous 
polynomial Q'^^{yi, y) of degree q and {ji^-v—q) in the components of the Jacobi coordinate 
vectors x and y, respectively. Q^''(x, y) is a solution of the Laplace equation. We have 
proved that — z/+ 1) generalized harmonic polynomials (5g''(x, y) constitute a complete 
set of independent base functions for the total orbital angular momentum [//, v] . Any 
wave function with the given angular momentum in the system can be expanded with 
respect to the base functions, where the coefficients are the functions of the internal 
variables, called the generalized radial functions. The three-body Schrodinger equation 
with a spherically symmetric potential V inD dimensions reduces to the generalized radial 
equations satisfied by the generalized radial functions. Only three internal variables are 
involved in the functions and equations. The number of both the functions and the 
equations for the given angular momentum [/x, v] is finite and equal to {ii — v Only 
a finite number of partial angular momentum states arc involved in constructing the 
generalized harmonic polynomials Qg'^(x, y), and the contributions from the remaining 
partial angular momentum states have been incorporated into those from the generalized 
radial functions. 

The generalization of this method to a quantum A'^-body system in £)-dimensions is 
straightforward. The difficulty is how to obtain the unified forms for the generalized 
harmonic polynomials, because it needs D — 1 vectors to determine the body-fixed frame 
and there are A^ — 1 Jacobi coordinate vectors. The cases with N < D are very different 
to the cases with N > D. We will study this problem elsewhere. 
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